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of event-by-event distributions of conserved quantities and lattice thermodynamic predictions based on the 
grand canonical ensemble. In this talk we discuss how the QCD crossover temperature and the freezeout 
curve are extracted from the analysis of fluctuations. We report that one can also go further and locate 
the QCD critical point at ~ 2Tc. We also list the systematics which must be brought under control in 
future. 
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1. Introduction 

In the last years a major development occurred in the experimental study of the phase diagram of QCD. It 
was established that certain measures of the shape of the distribution of event-to-event fluctuations of conserved 
quantities [1] which were predicted through lattice QCD simulations [2] agreed well with experimental observations 
in relativistic heavy-ion collisions [3]. 
The Maclaurin series expansion of the pressure, 

(where t = T/Tc and z = fi/T, T is the temperature, fj, the quark chemical potential and Tc any measure of the 
location of the QCD cross over at z = 0) is now the method of choice to construct various extrapolations to finite 
chemical potential on the lattice [4]. Eq. (1) is written in a form which emphasizes that dimensionless functions 
of dimensionless numbers are the output of lattice computations. The x^'^\t,z), which are the n-th derivatives 
of P with respect to ^ are called nonlinear susceptibilities (NLS); x*^^' is the quark number density and x'^' the 



* E-mail: sgupta&tifr.res.in 



Finite size scaling on the pliase diagram of QCD 



quark number susceptibility [5]. The series expansions of the NLS are obtained from eq. (1) by merely taking 
derivatives with respect to ^. Analysis of the series expansion of x'^''(ii z)/T'^ led to current estimates [6] of the 
position of the critical end point of QCD 



ZB = 1.8±0.1 and tjs = 0.94 ± 0.01, (2) 

where the choice made in [2, 6] was to use for the scale of the temperature, Tc(P), at which the Polyakov 
loop susceptibility peaks. The observables for which agreement between lattice predictions and experiments are 
demonstrated are the ratios [1] 



X^ '{i,z)/T^ X' '(i, 2)/T^ 

A development of this kind immediately makes possible other ways of looking at old questions and asking new 
physics questions; outlining some of these is one purpose of this review. The other is an equally important task: 
critically examining the assumptions that went into the comparison with a view to making them quantitatively 
testable. These are the contents of the next two sections. 



2. Old and new questions 

The initial demonstration of the agreement between experimental observation and lattice predictions proceeded in 
the following way. Lattice simulations predict dimensionless ratios as functions of other dimensionless ratios. In 
order to make contact with experiment two inputs were needed. The first was the scale of the lattice computations, 
Tc{P). It can be taken from current lattice measurements [7] to be 

7;(P) = 176±3±4 MeV, (4) 

where the first error comes from finite temperature lattice computations and the second from errors in the setting 
of scale by matching lattice computations to zero temperature hadron properties. Then lattice computations 
of x'"^(i,0) had to be resummed to give the ratios mi{t, z) along the freezeout curve, {Tf cor- 
responding to heavy-ion collisions performed with center of mass energy This second input, the freezeout 
curve, was taken from a parametrization of the hadron resonance gas model (HRGM) made through a fit to data 
on yields of various particles in [8]. 

2.1. Two classic questions 

It was pointed out in [2] that if reasonable agreement between lattice predictions and data were observed, then 
by relaxing the above input conditions, these measurements could be used to directly determine either Tc or the 
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Figure 1. Freeze out parameters at \/S — 200 GeV estimated from the STAR collaboration's measurements of mo and ma. The 
inner error bars show the effect of statistical errors in data propagated to this parameter (in some cases they are smaller 
than the size of the symbols), and the outer bars show the effects of systematic errors in data. The gray horizontal 
band is the statistical error band of the HRGM fit [8]. 

freezeout curve, or both. These are new ways of putting together lattice predictions and experimental data in 
order to examine classic questions in the field. 

The new program was initiated in [9] where the freezeout conditions were taken as before and the scale Tc was left 
free to be determined by the lattice predictions and experimental data taken together. Varying the scale Tc{P) 
in order to maximize the agreement between lattice predictions and data gave 

n{P) = 17512. MeV. (5) 

The errors are statistical only. The scale determination using the agreement of lattice predictions with either 
single hadron properties (eq. 4) or bulk matter (eq. 5, above) agree well, providing the first hard evidence that in 
lattice gauge theory we have a full theory of non-perturbative QCD. 

The second piece of physics that one can extract, following the program of [2], is the freezeout point at different 
x/S. In an expanding system net yields and fluctuations may freeze out at different points in the phase diagram (as 
we discuss in the next section), so this is an important measure of the degree of thermalization. The appropriate 
tools for this follow from eqs. (1) and (3), and some of the steps are given in [2]. These lead to the expressions 
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Here ri is the leading order estimate of the radius of convergence of the expansion of x'^V^^^i ^^'^ ''2 is the 
estimate at the next order. The t dependence of the right-hand side is hidden in ri_2 whereas the z dependence is 
explicit. Since the leading dependence on z is the factor of 1/z, these measurements have been called baryometers 
[10]. This is, of course, only approximately true (to within 5% at \^ = 200 GeV and to 25% at \/S = 62.4 GeV, 
as it turns out). In the following we will term the leading order (LO) approximation that is obtained by neglecting 
all the terms except 1/z. The next-to-leading order (NLO) approximation comes from keeping the second term 
in the series. Also shown is the resummation of these into a simple pole (the reason for doing this is given in [2]), 
an approximation that we term the NLOPade. 
Using the data from [3], the LO analysis yields 



= 200 GeV) = < 



0.1884 ± 0.0007 (±0.033) from mo, 

(7) 

0.183 ± 0.036 (±0.076) from mg. 



These estimates are shown in Figure 1. The first set of errors is statistical, and obtained neglecting covariance 
between measurements. The error estimates shown within brackets are similar, but obtained by adding statistical 
and systematic errors in quadrature. The error analysis should be treated as indicative, since it can be easily 
improved by experimental collaborations. The freezeout point inferred from HRGM fits is z — 0.142 ± 0.005, 
compatible with the above analysis when systematic errors are accounted for. Interestingly, the agreement between 
freezeout conditions from HRGM and fluctuations improves quite a bit in going from LO to NLO at \/S = 62.4 
GeV. 

Since there is no pure baryometer, for the NLO analysis one needs to specify t, or equivalently the values of ri_2. 
Given the near agreement of the freezeout value of z from yield and fluctuations analyses, it is would seem that 
no significant error is introduced if we take t = 0.94 as given in [8], instead of determining it through the longer 
procedure outlined in [2]. With this input and the lattice determination of [6] that ri = r2 — 1.8, one can perform 
the analysis to orders NLO and NLOPade. The results are compatible with the estimates in eq. (7) as shown in 
Figure 1. 

Clearly, there are systematic errors in the measurement of yields. However these are not taken into account 
in the extraction of freezeout conditions. It would be interesting in future to know how large these are, and 
their influence on the comparison on the two ways of determining the freezeout conditions. Furthermore, when 
sufficiently large number of cumulants of event-to-event distributions becomes available it may be possible to fit 
Tc and the freezeout conditions simultaneously, thus removing the necessity of external inputs. 

2.2. A new question 

The beam energy scan at RHIC is examining the exciting new question of the location of the critical point of 
QCD. The comparisons above show that systems produced at large \/S are thermalized. It is expected that if the 
systems are close to a critical point, then, due to increased relaxation times and correlation lengths, they will not 
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Figure 2. Contours of constant in the plane of mo and ma. The shaded area below the diagonal is not allowed for any real 
value of the critical point ze- The open curves above the diagonal are contours of ze for the values shown near the 
lines. The ellipses are Ict contours of the measurements at ^/S — 200 GcV and \/S — 62.4 GeV [3]. The outer ellipses 
are obtained when statistical and systematic errors are added in quadrature, and the inner ellipses show statistical 
errors only. 

be in local thermal equilibrium [11]. So the signal for the critical point will be that of agreement with thermal 
QCD away from the critical point and lack of agreement near the critical point. 

Interestingly enough, it is possible to confirm these results using only data far from a critical point by utilizing 
an interesting coincidence — that the freezeout value t{y/S — 200 GeV) — 0.94 is coincidentally very close to ts 
of eq. (2). At this value of t, the series expansion of mo and ms contain ri = r2 = ze- We have already used this 
in the above extraction of the freezeout value of z beyond LO. 

As a result of this coincidence, we can insert ri = r2 — ze in eq. (6), and then extract both z and ze directly 
from data. This requires a simultaneous fit of data on mo and ma at large -v/S", either top RHIC energy or LHC 
energies, to the NLOPade approximation. It is remarkable that under fairly weak assumptions one can extract the 
location of the critical end point directly from data at the highest collider energies without the direct intervention 
of lattice predictions. 

Using the same data which led to the fits shown in Figure 1 a simultaneous extraction of z and ze gives interesting 
results. Since the errors on mi are smaller, the best-fit value of z is close to that obtained in the NLOPade fit of 
mo. Simultaneous use of mo and ma also restricts the errors on z to be smaller than the maximum errors shown 
in Figure 1. 

As for the position of the critical point, in Figure 2 we show contours of constant ze in the plane of mo and ma, 
obtained by solving eq. (6). Whenever mo = ma one can invoke a Poisson description of the data and thereby 
push the critical point away to infinity. This special case is contained in the contour plots. Since almost all the 
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error ellipses cross the diagonal line in Figure 2, this crude analysis of errors cannot put an upper bound on ze- 
With this treatment we can only give very rough limits on the location of the critical point — 

ZE = 2.03i^°44 so that fiE = SbStTso MeV. (8) 

More precisely, one can bound ze > 0.59 at the 68% confidence level. This is in agreement with the lattice results 
quoted in eq. (2). 

Since it is clear from Figure 2 that the errors are dominated by those in ms, one can try to make estimates which 
take into account only the measurement of mo. Using as inputs the measurement of mo at VS ^ 62.4 GeV [3] 
and the freezeout z obtained from the fit of the HRGM [8], one gets 

ZE = 2.041°:;^ (1^29) , (9) 

where the first set of errors are statistical and the second from the composition of statistical and systematic. It 
is interesting that in either approach the best fit value of ze are similar, consistent with the lattice results at 2a. 
One expects improved results from the more careful error analysis that full access to data can give for two 
reasons — first, the obvious one that the since the publication of [3] much improvement has occurred in statistical 
and systematic analysis of data, resulting in decreasing error bars significantly; second, one expects a positive 
covariance of errors in mo and mg, especially the systematic errors. The net result is that the ellipse will shrink 
and tilt to the right, thus possibly giving good lower and upper bounds for ze- This presents a strong case for 
such analyses from RHIC and LHC experimental collaborations. 

3. Systematic errors and length scales in the fireball 

Systematic errors which afi'ect any comparison of theory and experiment can lie with either. The sources of 
systematic errors in lattice computations are well understood. We show there that there is some control over 
these, and the remainder can be brought under control with more computation. However, since the subject is still 
relatively young, understanding the experimental systematics involves understanding and exploring new physics. 

3.1. Lattice systematics 

If the theory were an ad-hoc model, one would not ask too much of it. However, when it is a first-principles 
method of computation in field theory one does need to examine systematic effects in various ways. Even in 
the comparison of perturbative QCD with collider data such questions are sometimes asked and occasionally 
the answers show the need for further work. For lattice thermodynamics the sources of systematic errors were 
enumerated in [12]; the main as-yet unquantified effect is of removing the lattice cutoff, i.e., the rate of approach 
to the continuum limit. 
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Figure 3. Lattice computations of the ratio 7713 [2] with two different cutoffs {Nt — 4 and 6). The two computations agree well 
with each other at \/S of about 50 GeV and higher. Below this the disagreement is due to cutoff effects in the location 
of the critical point. 



An exploration of this systematic effect [2] is shown in Figure 3. The result of a computation of ma with two 
different lattice spacings is shown. The predictions for \/5 > 50 GeV are clearly within control. There is less 
control at present at smaller lattice spacings. This can be traced to remaining lattice cutoff effects on the position 
of the critical point, which can be improved with further work on the lattice. This depends now largely on the 
computational power one can bring to the problem. 

3.2. Experimental systematics 

The most obvious systematic errors one must control in experiments is due to the fact that one is trying to apply 
the results of grand canonical thermodynamics to the data. One must then make sure that the experimental 
situation is close enough to local thermodynamic equilibrium: that the system is close to chemical equilibrium, 
that diffusive and advective phenomena are in balance, and that the part of the system under observation is 
neither too small nor too large. We will examine these questions in turn. 

Geometric acceptance cuts in experiments, such as those on centrality, the pseudorapidity and azimuthal angle, 
serve to define the volume of the fireball which is observed, Vot- If this volume is comparable to the full volume 
of the fireball, Vfb, then conserved quantities carry only initial state information. Fluctuations can still be 
observed, but they tell us about initial state fluctuations. These studies are of great interest because they can be 
correlated with other probes of the initial state [13]. However, none of the dedicated heavy-ion experiments have 
the acceptance to be able to do this. On the other hand the LHC detectors ALEPH and CMS with their near 47r 
coverage can easily be used to study initial state fluctuations of conserved quantities. 

It is necessary to have Vot ^ Vfb in order for the unobserved part of the fireball to act as a reservoir of energy 
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and particles so that thermodynamics in the grand canonical ensemble can be applied. Systematic corrections in 
powers of Vob/Vfb can be worked out [14]. Another necessary condition is that the scale of observation should be 
much larger than the microscopic length scale in the fireball. For fluctuations of baryon number, the appropriate 
microscopic length scale to examine is the longest static correlation length of the baryon quantum number. The 
inverse of this correlation length is called the nucleon screening mass, A4n- So, the condition that is usually 
required is Al^rVoi, 2> 1. At the critical point, where the divergence of x'^' implies the divergence of A4n, the 
inequality does not hold, and QCD thermodynamics is not expected to describe the fluctuations [f 1]. 
When the inequlity is valid, the observed volume contains many independently fluctuating volumes. As a result 
the shape of the distribution of fluctuations tends to a Gaussian; this is an application of the central limit theorem 
of statistics. The classic theory of fluctuations [15] relate the mean and variance of the distribution to quantities 
such as x'^' and x(2)- Systematic corrections in powers of A4%Vob relate higher cumulants to the NLS and lead 
to the study of observables such as those in eq. (3). The analysis of experimental data which was described in 
earlier sections is therefore part of the systematic theory of finite volume effects. 

There is another length scale which needs to be investigated. If the total number of baryons and anti-baryons 
seen in a given event is — B + B, then the volume per detected baryon in that event is Vob/B+ = i^'^. The 
distribution of ( is of some interest in understanding thermalization. Extreme values of are rare events in a 
thermal ensemble. If these appear more frequently than in a thermal ensemble, then that implies that initial 
random fluctuations in baryon number have not had time to equilibrate by diffusion. Such events would strongly 
influence tails of the event-to-event distribution of the net baryon number, and can distort measurements of the 
higher cumulants. Therefore, a study of high-order cumulants and their comparison with lattice predictions is of 
interest in understanding the speed of approach to equilibrium. It would be of special interest to divide the data 
sample into bins of(^<lor(;^>lor("^l, and compare the cumulants in different bins. 

The competing dynamics of diffusion and flow create other length scales in the plasma [16]. If the expansion rate is 
small enough then fluctuations are evened out by diffusion. However, if the expansion rate is larger, then at some 
scale fluctuations are frozen into the fiuid. Which of these actually occurs is discriminated by a dimensionless 
number called Peclet's number — 

£v 

Pe=|, (10) 

where £ is the length scale of interest, v is the flow velocity and D is the diffusion constant. When Pe 1 diffusion 
dominates and for Pe 2> 1 baryon number is passively transported by the flow. As a result fluctuations freeze out 
when Pe ~ 1. Using the fact that V ~ Cs/Mn, where Cs is the speed of sound, we find that Pe = M£Mn, where 
M is the Mach number. The freezeout scale for fluctuations, Ifo, is therefore 

iloMN ^ ^. (11) 



If the observations see near-equilibrium fluctuations, then clearly one must have Ifo < y'Vob- This implies that 
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Since there is no evidence for shock waves in the fireball, one can assume that AI < I, and 
therefore the observed volumes are sufficiently large for the usual finite size scaling theory to work. 
Chemical freezeout occurs when various reaction rates become slower than the expansion rate. Most strong 
interaction cross sections are of very similar magnitude, so treating them as having a common freezeout parameter 
is expected to be a reasonable approximation. The fact that a model, such as the HRGM, which is based on this 
assumption works as a good description of the fireball is then not surprising. 

ft was recently pointed out [17] that the isospin changing reaction, p -f-)- n which goes through the nearly resonant 
A channel has a very small energy denominator and therefore is relevant even after the chemical freezeout observed 
through fitted yields. This is an interesting observation since it can change the n/p ratio from its thermal value. 
Since neutrons are uncharged and therefore unobservable, estimates of baryon fluctuations are based on the 
assumption that proton and neutron fluctuations are identical within errors. This was checked through event 
generators in [3]. It would be useful to check how strongly the proposal of [17] affects event-generator estimates 
for proton fluctuations. 

4. Summary 

The observed agreement [3] of lattice predictions and experimental observations of fluctuations have many conse- 
quences. Among them are the possibilities of extracting measures of Tc or freezeout conditions from a comparison 
of data and thermodynamic predictions from lattice computations [2]. Interestingly, it also seems possible to ex- 
tract the location of the critical point indirectly from an analysis of fluctuations observed at the highest energies 
in RHIC and LHC (see section 2.2 and Figure 2). 

There are indications of good control over systematic effects in lattice computations (see Figure 3 and the sur- 
rounding discussion), and future work will settle this. Some systematic errors on experiments have been dealt 
with in detail before; others are discussed here. The fiuctuation observations of eq. (3) are finite size scaling 
quantities which take into account the effects of finite Vob-M'N- Methods which take into account non- vanishing 
Vob/Vfb have been discussed recently [14]. When Vot/Vft — 1 then one has a "persistence of memory effect" 
which allows the study of initial state fiuctuations in baryon number. The effect of fiuctuations in volume per 
baryon may affect high moments, but are amenable to experimental investigation. The competing effects of flow 
and diffusion give rise to a "Peclet length scale" [16] which controls the time at which fluctuations at the scale of 
Vob freeze out. Hydrodynamic computations can give more details of this process in future. Slow isospin changing 
reactions [17] have been noticed recently and their effects need to be understood. 

In summary, the study of fluctuations has become very interesting as it emerges that there is more to understand 
and much to gain by doing so. 
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